DESY 01-117 
WUE-ITP-2001-032 

Influence of CP and CPT on production 
and decay of Dirac and Majorana fermions 

G. MOORTGAT-PlCK a [], H. FRAAsf] 

a DESY, Deutsches Elektronen-Synchrotron, D-22603 Hamburg, Germany 
b Institut fur Theoretische Physik und Astrophysik, Universitat Wiirzburg, D-97074 Wiirzburg, 
Germany 

Abstract 

The consequences of CP and CPT invariance for production and 
subsequent decay of Dirac and Majorana fermions in polarized fermion- 
antifcrmion annihilation are analytically studied. We derive general 
symmetry relations for the production spin density matrix and for the 
three-particle decay matrices and obtain constraints for the polariza- 
tion and the spin-spin correlations of Dirac and Majorana fermions. 
We prove that only for Majorana fermions the energy and opening an- 
gle distribution factorizes exactly into contributions from production 
and decay if CP is conserved. 
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1 Introduction 



It is generally accepted that supersymmetry (SUSY) is one of the most 
promising concepts for physics beyond the Standard Model (SM). A special 
feature of SUSY models is the Majorana character of the neutralinos, the 
fermionic superpartners of the neutral gauge and Higgs bosons. After the 
observation of candidates for neutralinos at a future e + e~ collider their 
identification as Majorana particles is indispensable. Therefore an extensive 
study of the general characteristics of Majorana fermions produced in e + e~ 
annihilation is of particular interest. 

In ||, HI , the authors analyzed the consequences of CP and CPT imparl- 
ance for the symmetry properties of the cross section for Majorana fermions 
produced in e + e~ annihilation with polarized beams. They proposed useful 
methods to distinguish between Majorana and neutral Dirac particles via 
the energy distributions of the leptons from their leptonic decay. 

In the present paper we extend the investigation of CP and CPT sym- 
metry properties to the complete spin production density matrix for Dirac 
and Majorana fermions and to their decay matrices for a three-body decay. 
We consider the most general dependence on beam polarization so that our 
analysis is also applicable to annihilation with the exchange of Higgs 

bosons. 

CP and CPT symmetry relations lead to important consequences for 
the factorization of differential cross sections for production and subsequent 
decay into a production and a decay piece. In || it has been proven that 
the differential cross section factorizes only if the kinematic variables are 
properly chosen. However, factorization is ruled out in particular for the 
energy and angular distributions of the decay products in the lab frame due 
to interference between the helicity amplitudes. 

In numerical analyses || it was demonstrated that for production and 
three-particle decay of neutralinos in e + e~ annihilation the energy distri- 
bution of the decay products in the lab frame as well as the distribution 
of the opening angle between two of them indeed factorize. Contrary, for 
production and three-particle decay of charginos || the spin correlations 
between production and decay considerably contribute to the energy distri- 
bution of the decay products and therefore prevent factorization. Here we 
prove that the factorization of these observables in the case of Majorana 
fermions is a consequence of their CP/CPT symmetry properties. We show 
that the analysis of both the energy spectrum of the decay fermions and the 
opening angle distribution is useful for the identification of the Majorana 
character of neutral particles produced in fermion-antifermion annihilation 
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with polarized beams. 

This paper is organized as follows: In section 2.1 we give details of 
the spin density matrix formalism which is applied to the analysis of the 
consequences of CPT and CP invariance for the production density matrix 
of two different Dirac and Majorana fermions in section 2.2 and for three- 
body decay matrices in section 2.3. In section 3 we study analytically the 
consequences for angular distributions and the energy spectra of the decay 
fermions. 



2 Constraints on the production of Dirac 

and Majorana fermions from CPT and CP 

2.1 Spin-density matrix and cross section 

We consider pair production of Dirac or Majorana fermions in fermion- 
antifermion annihilation and their subsequent three-particle decay. The 
helicity amplitudes for the production processes 

/(pi, Ai)/(p 2 ,A 2 ) -> fi(pi,Xi)fj(pj,Xj) (1) 

are denoted by Tp^^ . In the case of Majorana particles it is fj = fj. For 
the decay processes 

fiipiAi) fa{m)fa{Pi2)fa{m), (2) 

fj{pj,Xj) -> fjl{Pjl)fj2{Pj2)fj3(m) (3) 

the helicity amplitudes are given by Tb^. and Here fa (fji) is a 

charged or neutral Dirac or a Majorana fermion and fa (fj2), fa (fj3) are 
Dirac fermions. The notation of the decays (||), (^) does not distinguish 
between fermions and antifermions for the initial particles and their decay 
products. The helicities of the decay products are also suppressed. 
The spin density matrices of the polarized beams can be written as 



P (f) = -(i+py) (4) 



1 

2 l 

p(i) = l^+py), (5) 



where Pj, Pj, Pj (Pj, Pp Pj) is the transverse polarization of / (/) in the 
production plane, the polarization which is perpendicular to the production 
plane and the longitudinal polarization. 
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The amplitude squared of the combined process of production and decay 
reads 



Ai A t A ■ A - 

1 J * J 



|T| 2 = £ IWtfWjtfpP^PnXtfDXv W 



It is composed from the (unnormalized) spin density production matrix 

Pp 3 ' XiXj = E pU)\\\Af)\' 2 \2 T PM\ 2 T ^{l^ ^ 

Ai A2A 1 A 2 

the decay matrices 

Pd,x% = Td ^ T dX ( 8 ) 

^D.A'.A, = T D,\ 3 T* DX , (9) 

and the propagators 

A(/ fc ) = l/\pl - ml + im k T k ], k = (10) 

Here pi, and denote the four-momentum squared, mass and total 
width of the fermion For these propagators we use the narrow-width 
approximation . 

We introduce for f, L (/j) three spacelike polarization vectors s a/i (/i) 
(s btJ, (fj)) which together with p^/mi (pj/rrij) form an orthonormal set. 
Then the spin density matrix of production and the decay matrices can 
be expanded in terms of Pauli matrices a a with the first (second) row and 
column corresponding to the helicity A = | (— |) M: 

Aj At , A - A 



pp 3 " 3 = ^A;^A'.nM) + 5 AjA 'E^A; s p^) 

3 a=l 

+^a : E<a'^) + E ^<^ a p(f^ en) 

b=l J a,6=l J 

^ = ^a^(/0 + E< :a S^(/,), (12) 

a=l 



6=1 



Here a (6) = 1, 2, 3 refers to the polarization vectors of fi (fj). The contribu- 
tions ££(/;) (£^(/,)), £&(/<) (£$>(/,■)) are linear, and is bilinear 
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in the polarization vectors s afJ- (fi) (s bu (fj)). In (11) the dependence of pp 
on beam polarization has been suppressed. 

The polarization vectors s 1 , s 2 , s 3 form an orthogonal right-handed 
system in the lab system: 

• 58 (fi) ( s3 (fj)) is i n the direction of momentum pi (pj), 

• s2 (fi) = |^^' = s 2 {fj) is perpendicular to the production plane, 

• s l {fi) (s 1 (/j)) is in the production plane orthogonal to the momentum 
Pi (Pj)- 

Then Ep(/y)/P(/j/j) is the longitudinal polarization, ^ P (fi,j)/P(fifj) 
is the transverse polarization in the scattering plane and %p(fi,j) / P{fifj) is 
the polarization perpendicular to the scattering plane. The terms Sp (fifj) 
are due to correlations between the polarizations of both produced particles. 

The amplitude squared |T| 2 of the combined process of production and 
decay (|6|) can be written as: 

iti 2 = ^A^piA^oi^m^o^^^^o + Espaos^^)^^) 

a=l 

+j2^ b p(fj)^D(fj)D(fi) + E xfifiMzumUfj)], (14) 

6=1 o,6=l 
and the differential cross section is 

da = -L|T| 2 (2vr) 4 ( 5 4 (p 1 + p 2 - E^hps, (15) 

AS 

I 

with the Lorentz invariant phase space element (flips. 

For the case of neutralinos and charginos the complete analytical expres- 
sions for the production density matrix and for the decay matrices are given 
in different presentations in || ||]. 

2.2 CPT and CP symmetries of the spin— density matrix 

In this section we derive constraints from CPT and CP invariance for the 
density matrix ( JTl] ) for the production of two (charged or neutral) Dirac 
fermions and of two different Majorana fermions, respectively, for polarized 
beams. 

All contributions from the exchange of particles a, (5 to the coefficients 
£ = {P,S P ,Sp,S P 6 } (16) 
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of the spin-density matrix ( JTT| ) are composed of products C of couplings, 
the propagators A(a), A(/3) and complex functions S of momenta and po- 
larization vectors: 

£ ~ i?e{C x A(a)A((3)* x 5}. (17) 

Assuming complex couplings and taking into account the finite widths of 
the exchanged particles the general structure is 

£ ~ Re(C)Re[A(a)A((3)*]Re(S) - Re(C)Im[A(a)A((3)*]Im(S) 

+Im(C)Re[A(a)A(P)*]Im(S) - Im(C)Im[A(a)A(py]Re(S), (18) 
= £rrr + £rii + £iri + £iir- (19) 

If CP is conserved all couplings can be chosen real ]9j and £jri = £iir = 0. 
The terms £ru and £hr are due to interference between s-channel exchange 
of particles with finite width and the crossed channels. Their contributions 
can be neglected except in the vicinity of the s-channel resonances (exchange 
of gauge bosons in the case of e + e~ or qq annihilation and in addition of 
Higgs bosons in fJ> + n~ annihilation), since they are proportional to the width 
of the exchanged particles far from the pole. In £ru and £iri the imaginary 
part of S originates from products of momenta and polarization vectors 
with the Levy-Civita tensor e^ upa - These two terms lead to triple product 
correlations of momenta which are sensitive to CP violation || • 

2.2.1 CPT invariance 

Under CPT the helicity states of a Dirac fermion fo and a Majorana fermion 
/m, respectively, transform as 

\f D (p,X) > C 4 T (_i)A-!|/ ( p -;-A) >, (20) 

|/m(p, A) >°^ T (-1) A -^ CPT |/m(p, -A) > , (21) 

with the CPT phase r] T = ±i of Majorana fermions |Kj. Beyond that 
time reversal would interchange in- and outgoing states. However, if the 
finite width of the exchanged particles is neglected CPT invariance and 
the unitarity of the S-matrix leads in leading order perturbation theory 
to symmetry relations for the helicity amplitudes and for the production 
density matrix. 

In the following we derive these symmetry relations from CPT invariance 
for (charged or neutral) Dirac fermions and for the special case of Majorana 
fermions. 
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2.2.1 



1 Dirac fermions 



CPT invariance relates in lowest order perturbation theory the spin density 

A j Aj , A' A' 

matrix p p 1 3 of the process 

// •/,(//,. A, A,) (22) 

A i A j , A' A' 

to the spin density matrix p p 1 1 of 

// •/,(//,. A, A, ! (23) 

by 

/V A " A ' A ' (P, 1 ' 2 ' 3 ; P}' 2 ' 3 ; 9) = „,/ V A • a! {p}\ -P|; P}\ -P|; 0)*, (24) 

where denotes the angle between / and f\ in ( |22| ) and = 7r — is the 
angle between the direction of beam particle / and the outgoing antifermion 
fi of the conjugated process 
larizations of / and / are interchanged with sign reversal of the longitudinal 
polarization. To derive ( p4|) the CPT transformation has been supplemented 
by a rotation Tv!^ 71 ") around the normal to the production plane so that the 
beam direction is unchanged (see Fig. |l]a). 

Both density matrices pp and pp can be expanded according to ( |Tl"| ) 
with the coefficients P(fifj), £?>(/>), & P (fj), ^f(fifj) and P(fifj), 
^p(/j)> ^-''p(fifj)- Then from CPT and ^(tt) invariance and with the 
substitutions 



(p l,2,3. p l,2,3. @) ^ (p l, 2j _ p 3. p l, 2) _p3. q = n _ @)) (25) 

one obtains the following relations: 

P(fif j ) = P(fif j ), (26) 

= r/aS^C/i), (27) 

^ b p(fj) = ribZpifj), (28) 

Xf(fif j )=VaVbV P b (fif j ), (29) 



with 771 = 1 = r/2 and % = — 1. 



23) in the CMS. For polarized beams the po- 



7 



2.2.1.2 Majorana fermions 



Since it is pp = pp in (24) the invariance under CPTx^l 71 ") leads in lowest 



order perturbation theory to constraints for the production density matrix. 
The different terms in the expansion ( |TTj ) of pp can be classified into con- 
tributions which are symmetric (Type S) and antisymmetric (Type A) with 
regard to the substitution 

(p l,2,3. p l,2,3. @) ^ (p l,2 ; _ p 3. pl^ _pj. ir _ @)) (30) 

where denotes the angle between / and fi in the CMS. 

• Type Si P \ ^p^ ^p^ ^p ? ? 5j*p^ . 5 ^.p 

• Type A: T,p, T.f, 

All of them include terms £rrr and £irj of ( ^8| ) . 

For unpolarized beams the terms of type S are forward-backward sym- 
metric whereas those of type A are forward-backward antisymmetric. These 
symmetry properties hold also for production of Majorana fermions with po- 
larized e + e~ beams. 

The forward-backward symmetry of the term P results in the FB- 
symmetry of the differential cross sections for the production of Majorana 
fermions with unpolarized and longitudinally polarized beams @. Beyond 
that, also the symmetry properties of the different components of their po- 
larization and of the spin-spin correlations are specific for their Majorana 
character. 

2.2.2 CP invariance 

We now study the consequences of the invariance under a CP transformation 
(illustrated in Fig. |l|b) for the production of two Dirac or two Majorana 
fermions. 

Under CP the helicity states of a Dirac and a Majorana fermion trans- 
form as 

|/ D (p,A)>^|/ D (-p,-A)>, (31) 
\f M (p, A) >^ rj CP \f M (-p, -A) > (32) 
with r/ CP = ±i. 



S 



2.2.2.1 Dirac fermions 

CP invariance results in a relation between the spin density matrix pp for 
the conjugated process 

// -> fi{Pi,h) fj (.Pj Aj) (33) 
and the density matrix pp for 

ff^fM,\)fM^j) ( 34 ) 



pT"^ i/'r- *f*\ 8) = P r V A • A A ' (P/, -Pf ; /'/• /f: 6), (35) 



with the same denotation as in d24|). For polarized beams the polarizations of 
/ and / are interchanged with an additional sign reversal for the longitudinal 
polarization and the polarization perpendicular to the production plane. 

Expanding both density matrices pp and pp according to (|ll|) one ob- 
tains with the same substitutions 

(P 1 / 2 ' 3 ; Pj' 2 ' 3 ; 9) -> (Pj, -P|' 3 ; P), -P 2 ' 3 - 9 = vr - 9) (36) 

as in (35) the following relations between the coefficients: 

P(fif j ) = P(fif j ), (37) 

= ^(/i), (38) 

E^(/ i )=r /6 E|,(/,), (39) 

E a P b (fJ j )= Va r ]b E a p b (f i f J ), (40) 

with 7/i = +1 and r/2 = — 1 = 7/3. For unpolarized or longitudinally polarized 
e + e~ beams, the dependence on the beam polarization is given by the two 
factors (1 — PfPj) and (Pj — P|)- If the widths of the exchanged particles 
in the s-channel are neglected, the CPT symmetry relations (|24]) hold and 
one obtains 

Ep = 0, Ep 2 = = Ep 1 , Ep 3 = = Ep 2 . (41) 



2.2.2.2 Majorana fermions 



With pp = pp in fl35|) CP invariance again leads to a classification of 
the different contributions in the expansion (|ll]) of pp according to their CP 
symmetry properties with regard to the substitution 

( pl,2,3. p l,2,3. @) ^ _ p 2,3. p ^ _p2,3. ^ _ ^ (42) 

where 9 is the scattering angle in the CMS: 
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• Type S: P, Sp, Sp 1 , Ep 2 , Ep 3 , Ep 3 , Ep 2 , 

• Type A: Sp, Ep, Ep 2 , Ep 1 , Ep 3 , Ep 1 

All terms have contributions with the structure £rrr and £rh of (|l8|). 

For unpolarized and longitudinally polarized e + e~ beams, the terms of 
type S are FB-symmetric and those of type A are FB-antisymmetric. For 
both Dirac and Majorana fermions the terms Ep, Ep, Ep 1 , Ep and Ep 2 
vanish for unpolarized or longitudinally polarized beams if the CPT relations 
(p6|)-(ps|) hold. The FB-symmetry or antisymmetry of the polarization and 
of the spin-spin correlations, however, is specific for the Majorana character 
of the produced fermions beyond the FB-symmetry of their production cross 
section. Since in general the production of Dirac fermions will not exhibit 
these symmetry properties the measurement of polarization and spin-spin 
correlations are helpful for an experimental verification of the Majorana 
character. 



2.2.3 CP violation 

If CP is violated due to complex couplings the two types £iri, £iir of CP 
violating terms contribute in ( |l8| ) , whereas far from the s-channel resonances 
£hr can be neglected as explained in section [2.2|. 



Since for Dirac fermions the CP transformation relates the two processes 
// — > fifj and // — > fifj we restrict ourselves to the discussion of the 
consequences of CP violation for the production density matrix of Majorana 
fermions for unpolarized or longitudinally polarized beams. Both types of 
CP violating terms £iri and £hr lead to contributions which show different 
FB angular dependence compared to the CP conserving ones. 

In particular terms with the structure £iri result in nonvanishing polar- 
ization Ep perpendicular to the production plane and nonvanishing spin- 
spin correlations Zj p and h p : 

• FB-symmetric in £iri: Ep, Ep, Ep 1 

• FB-antisymmetric in £iri: Ep , Ep 2 

These terms are proportional to Im(S) and are sensitive to triple product 



correlations of momenta. They are CP-odd quantities [11]. 



All other contributions in (11) are CP-even quantities and get in the 



case of CP violation terms with the structure £ur: 
• FB-symmetric in £ur: Ep, Ep 3 , Ep 1 
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• FB-antisymmetric in £ IIR : P, £p, Ep 1 , Ep 2 , Ep 3 

The consequences of CP and CPT for the FB-symmetry of the different 
terms of the production spin density matrix of Majorana fermions are sum- 
marized in Table |l|. 

2.3 CPT and CP symmetries of the decay matrices 

In this section we derive from CPT and CP invariance constraints for the de- 
cay matrices for the three-body decay of neutral and charged Dirac fermions 
and of Majorana fermions. 

2.3.1 CPT invariance 



2.3.1.1 Dirac fermions 

For the three-body decay 

fiifi, \) -> fil(Pil)fi2(Pi2)fi3(m) (43) 

of a Dirac fermion /j into a Dirac fermion fn and a fermion-antifermion pair 
fi2fi3, CPT invariance relates the decay matrix pp of (f43|) with the decay 
matrix po of the process 

fiipi, A;) -> fii{pix)fi2{m)fa{m)- (44) 

Since we do not study the polarization of the decay fermions we sum about 
their helicities and obtain 

with the momenta of the fermion-antifermion pair interchanged and the 
signs of the helicities reversed. 

Expanding pjj and pD according to (|l^) leads for the coefficients to 

D(m,m) = D(p i3 ,p i2 ), (46) 

x a D (m,m) = -t a D {m,m), a = 1,2,3. (47) 



2.3.1.2 Majorana fermions 



For the decay 

fi (p, \) -> filiPil) fi2 (Pi2 ) Ii2 (m ) (48) 
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of a Majorana fermion /j into a Majorana fermion fa and a Dirac fermion- 



antifermion pair fafa with = pd in (45), CPT invariance leads to 



constraints for the decay matrix from which one obtains definite symmetry 



properties for the expansion coefficients of pp in (12) 

D{pi2,m) = D(p i3 ,p i2 ) (49) 
^Um,m) = -Z a D {p l3 ,p i2 ), a = 1,2,3 (50) 

if the momenta of the fermion-antifermion pair fafa are exchanged. 
2.3.2 CP invariance 

Applying the space-inversion P to the decay ( |43| ) also the momentum of the 
decaying particle fa is reversed. We therefore add a rotation IZ^i:) around 
the normal to the plane (pi,pa). Then VxTZa(7r) is equivalent to a reflection 
at the plane (pi,pa)- 



2.3.2.1 Dirac fermions 



The combined transformation CP x IZ^tt) leads for the decay of a Dirac 
fermion to the process 

fiifi, ~A») -> fil{pil)fa{Pi2)fa(Pi3), (51) 

where pj2 and denotes the momentum jfo and pt 3 , respectively, reflected 
at the plane (pi,pn). The invariance under CP x TZft(^) relates the decay 
matrix pp for (|43|) with the decay matrix p£> for the antiparticle in (|5l|) 

P DA ' A fe,te) = (-l) Al ~ Al P D _ A '_ A .(Pi3,K2)- (52) 



In contrast to the CPT relation (45) the momenta of the fermion-antifermion 
pair are interchanged and additionally reflected at the plane (pi, pa). For 
the coefficients in the expansion ( |l2|) one obtains 

D{pi2,m) = D(p i3 ,p i2 ), (53) 

^D(m,m) = £z>(&3, Pi2), (54) 

S^ 3 (pi2,K3) = -S^ 3 (p i3 ,Pi2)- (55) 



Combining these CP relations with the CPT properties (|4q)-(p7D of the 
expansion coefficients leads to 

• symmetric D, S^ 3 

• antisymmetric Y? D 
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under the reflection of the momenta pa, and Pa at the plane (pi,pn)- 



2.3.2 .2 Majorana fermions 

For the decay (|48| ) of a Majorana fermion with pD = p^> in ( f32"|) the com- 



bined transformation CP x TZ^tt) constraints the decay matrix and their 
expansion coefficients when the momenta of the Dirac fermion fa and an- 
tifermion fa are interchanged and reflected at the plane of the Majorana 
fermions fi, fa. Neglecting the width of the exchanged particles CP and 
CPT invariance lead to the same symmetry properties of D, S^ 2 ' 3 as in the 
case of Dirac fermions under the reflection of the momenta pii, Pa at the 
plane (pi,pn). 



3 Energy and angular distributions of the decay 
products 

In this section we study for the production of Dirac and Majorana fermions 
and their subsequent three-particle decay into fermions the influence of the 
polarization and of the spin-spin correlations on the energy spectrum and 
angular distributions in the lab frame. If the decay of only one of the pro- 
duced fermions, e.g. fi, is considered, it is in (14) T> b D (fj) = and D(fj) = 1. 
Then the total cross section for the combined process of production and de- 
cay is given by a P (ff -> Ufj) x BR(fi -> fa fa fa). 



3.1 Energy spectrum and opening angle distribution of the 
decay leptons 

We expect that in general the energy distribution of the decay particles in 
the lab system depends on the polarization of the decaying particle. In its 
rest frame the decay angular distribution is determined by the polarization 
with respect to a suitably chosen quantization axis. Since boosting in this 
direction into the lab frame the energy of a decay particle depends on the 
orientation of its momentum with respect to this axis it also depends on 
the polarization of the decaying particle. The same argument applies to 
the distribution of the opening angle in the lab frame between the decay 
products. 

In the following we show that due to the specific CPT/CP properties 
of Majorana fermions the energy distribution and the opening angle dis- 
tribution of the decay fermions are independent of the polarizations of the 
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decaying particle. In this case both distributions factorize into a production 
and a decay piece. 

For our analysis we parametrize the phase space in the lab frame in a 
way that it factorizes into the production and decay phase space. 

3.1.1 Energy spectrum 

For the energy distribution da /dEi 2 the phase space can be parametrized by 
the polar angle 0^2 (@i,i3) between the momenta of fa and fai (fa) and the 
azimuth ($1.13) of fa (fa). All polar (azimuthal) angles are denoted 
by @ a f3 ($0/3) > where the first index denotes the polar axis: 

da = T\T\ 2 sin QdQdQ sin Bi^dQ^d^i^ sin Qi^d&i^d^i^dE^, (56) 

where is the production angle. For light fermions fa and fa$ T is given 
by: 

T = Q E i2 [mf - m\ - 2E i2 (Ei - gcos0 M2 )] ,^„, 

2H (2Tr)7 unTiE* [E, - </cos0 M3 - E i2 (l - cos 6^3 )] 2 ' 1 ' 

E), denotes the beam energy, E a denotes the energy of the fermion f a in the 
lab frame and q = \f>i\ = \pj\- The opening angle 0j2,i3 between fi 2 and fi3 



can be expressed via (A.l) by the independent variables chosen in (pq). 



The kinematical boundaries for 0^2 depend on Ei 2 : 

S 9 9 

[0,7r], o<e i2 <- 



©i,i2 = \ , , F m 2 -m 2 ^^ m 2 - m 2 ' "m 2 - m 2 (•*&) 

[ (0,arccos(^- j^,; 1 )). ^-+t<^ 2 <^-^. 
3.1.2 Opening angle distribution 

For the distribution of the opening angle 0i2,i3 between the fermions fa and 
fa from the decay fa — > faifa 2 fi3 it is favorable to parametrize the phase 
space of /j3 by 0j2,i3 and the corresponding azimuthal angle $12,13 (Fig. Q), 
i.e. to substitute in ( |56"| ) ©^3 — > 0i2,i3 and $^3 — > $12,13 ■ The kinematical 

9 9 

region of depends on i)i2 : < E i2 < 2(E i -qco S e i i2 ) ; T lie an S les ®M3 
and $1^3 can be expressed by ( |A.3|) and (|A.6|) , (|A.8| )-( A.1C| ) by the chosen 
independent angular variables. 
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3.1.3 Transverse polarizations 

In the contributions Sp(/j)S| ) (/j) ! a = 1,2,3, from the polarization of the 
decaying particle /j both factors depend on the polarization vectors s afJ, (fi) 
defined with respect to the production plane, cf. section 2.1. They can be 
written 

E£(/i) = and = E^»(/i). (59) 

To study their influence on the energy and the opening angle distribution we 
distinguish between the contributions Sp 2 (/j) from the transverse polariza- 
tions and S|>(/j) from the longitudinal polarization. In the contributions of 
the transverse polarizations, the coefficients E^ 2 (/j) depend on the azimuth 
&i,i2 between the production plane and the plane defined by the decaying 
fermion /j and the decay product fa. To separate the $^2 dependence we 
introduce a new system of polarization vectors t au (fi) with t 3l/ (fi) = s 3v (/j), 
whereas in the lab system t 2 (fa is perpendicular to the plane defined by the 
momenta of ft and fa and is in this plane orthogonal to the momen- 

tum of the decaying fermion f{. Then in ( |59"|) the transverse polarization 
vectors t >' 2v {fi) are 

s^ifa = cos^t^ifa-sm^t^ifa, (60) 
s 2u (fi) = sin^ot^C/O + cos^o^C/i). (61) 



Applying the transformations flP|) and (|Aj), flO|) , dATsj) - (]A.10Q for 



the angles 0^3 and $^3 to the phase space parameterization ( pq) the con- 
tributions of E^ 2 (/j) to the opening angle distribution and to the lepton 
energy spectrum vanish for both Majorana and Dirac fermions due to the 



integration over $^2 |12 |. 



3.1.4 Longitudinal polarization 



3.1.4 .1 Majorana fermions 

The longitudinal polarization £p(/j) is forward-backward antisymmetric if 
CP is conserved, cf. section [2.2.2| . Due to the factorization of the phase 
space in production and decay also the contribution of the longitudinal po- 
larization vanishes after integration over the production angle G. 

Consequently both the energy and the opening angle distribution of the 
decay products of Majorana fermions in the laboratory system are indepen- 
dent of spin correlations and factorize exactly into production and decay if 
CP is conserved. 
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The factorization of the energy and opening angle distribution is es- 
sentially equivalent since for both observables it is a consequence of the 
FB-antisymmetry of the longitudinal polarization of Majorana fermions. 

Assuming CP violation and taking into account the widths of the ex- 
changed particles the exact factorization of the energy distribution and of 
the opening angle distribution of the decay products of Majorana fermions is 
destroyed, since the longitudinal polarization Xp gets FB-symmetric contri- 
butions from CP violating terms. These additional terms have the structure 
Sim, and are therefore for energies far from the resonance proportional to 
the width of the exchanged particle. 



3.1.4.2 Dirac fermions 



In general, the longitudinal polarization of produced Dirac fermions is not 
forward-backward asymmetric. Therefore the energy spectra and opening 
angle distributions of the decay products do not factorize into production 
and decay. 



3.2 Decay lepton angular distribution 

For the distribution da/dcos Qi^2 of the angle Si,j2 in the lab frame between 
the incoming fermion / and the fermion /j2 from the decay fi — > fnfi-zfts 
we parameterize the phase space by 

da = F\T\ 2 sin <dd®d<$> sin Qi^dQi^i^ sin Qi^dQ^d^dE^. (62) 

and express all angles in T (|57| ) and the azimuth $^2 in (|60|), ( |6l"l) by the 
independent variables in (|62|). The relations are given in ( |A.2j ), QA.4 ), flA.5| ), 



|A.7| ) [12]. Thus the kinematic factor T depends explicitly on the scattering 
angle and neither the contributions (fi) of the transverse polarizations 
nor that of the longitudinal polarization Ep(/j) vanish due to phase space 
integration. 

Consequently neither for Dirac fermions nor for Majorana fermions the de- 
cay lepton angular distribution in the lab frame factorizes in production and 
decay but depends sensitively on the spin correlations. 



3.3 Siamese opening angle distribution 

The Siamese opening angle 0j2,i2 denotes the angle between decay products 
fi2 and fj2 from the decay of different particles fi and fj (fj in the case 
of Majorana fermions) in the lab frame. The distribution da/dcos@j 2 ^ 2 
is determined by the spin-spin correlations £p (|i~4] ) between the decaying 
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mother particles /, and Here it is favorable to parameterize both phase 
spaces for the decay of fi and of fj by the angles 0^3, $i,«3, and 0ij 2 , *&i,j2 

da = TQ\T^ sin QdQdQ sin Qi^dQi^d^i^ sin Qij3d&i t j3d^ij^dEj2 

x sin 0i,i3^0i,i3^i,i3 sin Qj2,i2d@j2,i2d^j2,i2dE i2 , (63) 

where 0ij3, $ij3 is defined with respect to the direction of fi and 0j2,i2 5 
^i2,i2- T is given by @) and 

1 £ j2 [m| - m|i - 2£ j2 (Ej - q cos 8j J2 )] 
2 5 (2vr) 5 m i r j [E, - q cos iJ3 - E j2 {l - cos j2 , i3 )] 2 ' 

Using the transformation formulae (A.l) and ( [A, 11 ) - ( |A.16[ ) all angles can 
be expressed by the independent variables chosen in (|63]), 

With the same arguments as for the contributions from transverse po- 
larization to the opening angle distribution one infers: 

The Siamese opening angle distribution factorizes neither for Majorana nor 
for Dirac fermions, it depends, however, only on the diagonal spin-spin cor- 
relations Ep 1 , Sp, Sp if CP is conserved. 

The polarizations of the produced fermions and the spin-spin correla- 
tions between these contributing to the energy and the different angular 
distributions of their decay products are listed in Table ^. 



4 Conclusion 

Assuming CPT and CP invariance we derived symmetry properties of the 
spin density matrix for production of Dirac and Majorana fermions in fermion- 
antifermion annihilation with polarized beams and for their three-particle 
decay matrices. Majorana fermions show specific forward-backward sym- 
metry properties of their polarizations and the spin-spin correlations. In 
particular for the production of Majorana fermions with unpolarized or lon- 
gitudinally polarized e + e~ beams their longitudinal polarization in the CMS 
is forward-backward antisymmetric so that the energy distributions of the 
decay products and the distribution of the opening angle between them are 
exactly independent of spin correlations if CP is conserved and if the width 
of the exchanged particles can be neglected. Thus they factorize into a 
production and a decay piece which allows to study the dynamics of the de- 
cay of Majorana fermions independently of that of the production process. 
Since this factorization is specific for Majorana fermions it opens the possi- 
bility to establish experimentally the Majorana character by comparing the 
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measured decay energy and opening angle distributions with Monte Carlo 
studies. 
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A Parametrization of phase space 

For completeness we list all angular relations used for the different parametriza- 
tions of the phase space in section 3. 

cos 6 i2 ,i3 = cos 8^2 cos 0^3 + sin 8 i)i2 sin Q iji3 cos($ iii2 - $i,i3), (A.l) 
= cos 81,12 cos 6 M3 + sin 8^2 sin 81^3 cos(3>i ii2 - $1,13), (A.2) 
cos6 i)i3 = cos 8 i)i2 cos 8i2,i3 + sin 8^28^812,13 cos $ i2 ,i3, (A. 3) 

= cos cos 81^3 + sin 6 sin 61^3 cos $1^3, (A. 4) 

sin©ii3 . 

sm$ii 3 = . ' sm$ia, (A. 5) 

sm Q i:i3 

= sin $ iji2 cos($ iji2 - $1,13) - cos $ iji2 sin($i >i2 - $^3), (A. 6) 
cos$ iii3 = -[cos 81,13 sin 8 - sin 81^3 cos 8 cos $i )i3 ]/ sin 8^3, (A.7) 
= cos $ M2 cos($ M2 - $ iii3 ) + sin sin($ M2 - $i,i3)- (A.8) 

In (|A.6|) and ([A .8]) one has to insert: 



, \ COS 0i 2 ,i3 Sm 8i i 2 - COS 8j i2 Sm 8 i2 ,i3 COS $ i2 ,i3 ,. n s 

cos($ M2 - $ i i3 ) = ! ; . 77 ■ — , (A. 9) 

sm i>i3 

sin($ M2 - $ M3 ) = sm ® i2 - i3 sm$ i2ti3 . (A.10) 

sm 8 M3 

The angles 81,12, $1,12 are given analogously to ( |A.4| ), (|A.5f) and ( |A.7| ) with 
8l,i3 -> ©1,12, *i,t3 -» $l,i2 and 0j )i3 -> 0j >i2 . 

COS 0^,2 = COS 0jj2 COS 0j2,i2 + SU1 0jj 2 sin 0j 2 ,i2 COS $j2,i2, (A. 11) 

sin$ i)i2 =- sm($ i)i2 - *»,»2) cos$ iJ2 + cos($ iij2 - $i, i2 ) sin $ iij2 , (A.12) 
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COS $^2 = Sin($jj 2 - ^ > i,i2)sin$j J -2 + COs($jj 2 - ^i,i2)cOS$jj2, (A.13) 

cos &j2,j3 = cos &i j 2 cos Qij 3 + sin 8^2 sin 8^3 cos($;j 2 - ®i,j3), (A. 14) 
where one has to insert: 

sin($ij 2 - *i,o) = Sm ®f '* 2 sin * i2 ,«, (A. 15) 

sin fe>i,i2 

cos($ij2 - *i,»2) = — . p — — ■ (A.16) 





cp, r a 


= 


Zjp — 


, p vi vii v 22 v^3 

o. z , ^p-i ^ P 5 P ' P 

q • v 3 V i3 v 3 i 

et. Zjp, ^ p i ^ p 

0,S P 2 = = S P 1 ,S 2 3 = = S P 2 


QP, T a = 


= 


cp,r Q ^o 


c . p vi v 2 vii v 22 

o. JT , Zjp , Zjp , .P ? jP 

p . v 3 V 13 V 31 

d. ^ P ) P 5 


y33 yl2 y21 
, Zjp , Zjp , Zjp 

y23 v32 
Zjp , Zjp 


„. p vl VH v 22 v 33 v 23 v 32 

o. J , Zjp, Zjp j Zjp j Zjp j Zjp , Zjp 

V 2 V 3 V 12 V 2i V i3 V 3i 

d. Zjp , Zjp, Zjp , Zjp , Zjp , Zjp 



Table 1: The forward-backward symmetry (s) and antisymmetry (a) of 
all terms of the production spin-density matrix (|ll|) of Major ana fermions 
for unpolarized or longitudinally polarized e + e~ beams with (r a ^ 0) and 
without (T a = 0) consideration of the width of the exchanged particles a 
for CP conservation (CP) or CP violation (C/P). 





Dirac fermions 


Major ana fermions 


Decay distrib. 


CP 


qp 


CP 


QP 


energy 


E P 


Zjp 




yA 
Zjp 


opening angle 


Sp 


Y 3 
Zjp 




y3 
Zjp 


angular 


vl,2,3 
Zjp 


v l,2,3 
Zjp 


v l,2,3 
Zjp 


v-,1,2,3 
Zjp 


Siamese angle 


v ll,22,33 




v ll,22,33 
Zjp 


\^ab 
Zjp 



Table 2: For the energy and the different angular dis- 
tributions of the decay fermions in the lab frame the 
contributing polarizations and spin-spin correlations 
are specified. For CP conservation (CP) the polariza- 
tion dependence is different for Majorana and Dirac 
fermions. If CP is violated (QP) all distributions de- 
pend on spin correlations. 
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/(pi.Ai) 





e 



f{P2,M) 



b) 

/(pi,Ai) 





/(P2, A 2 ) 



fi{pu 



CPT 



/(-P2, -Ai) 



CP 




/(Pi,-A 2 ) 



fi(-Pi,-Xi) 



Figure 1: Production processes in the CMS under a) CPT followed by a 
rotation 7^-2 (tt) and b) CP transformation. 




Figure 2: Definition of momenta and polar angles in the lab system. The 
indices of the angles denote the plane covered by the corresponding mo- 
menta, the first index denotes the corresponding polar axis. The momenta 
and angles in the decay of fj are denoted analogously. 
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